Fronts dynamics in the presence of spatio-temporal structured noises 
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Front dynamics modeled by a reaction-diffusion equation 
are studied under the influence of spatio-temporal structured 
noises. An effective deterministic model is analytical derived 
where the noise parameters, intensity, correlation time and 
correlation length appear explicitely. The different effects of 
these parameters are discussed for the Ginzburg-Landau and 
Schloegl models. We obtain an analytical expression for the 
front velocity as a function of the noise parameters. Numeri- 
cal simulations results are in a good agreement with the the- 
oretical predictions. 

PACS: 02.50.-r,05.40.-a 



I. INTRODUCTION 

The role of external fluctuations in extended systems 
is a subject of very active research because of its rele- 
vance in pattern formation in non-equilibrium systems 
JjJU ■ A simple example of a dynamical pattern is a front 
moving at constant velocity. Fronts can be easily mod- 
elized by a reaction-diffusion equation with two steady 
states of different stability ||. The study of the front 
dynamics under the influence of noises is relevant not 
only from theoretical point of view but also 

from practical point as recent works on chemical kinetics 
have shown |Tq,0,0,E3|. In these experiments a chemi- 
cal wave moves under the influence of and external fluc- 
tuating illumination which is projected in the reactive 
medium. This external source of noise has finite inten- 
sity, correlation time and correlation length. Thus, in this 
system we face with spatio-temporal structured noise and 
not with a white noise. The present work could be use- 
ful to clarify the role of these parameters on propagating 
structures as studied experimentally in |fiof . 

Previous studies dealt with this problem under the sim- 
plified assumption of white (<5-correlated) external fluctu- 
ations 

gUJZlffi. Nevertheless one can ask about the 
correctness of this assumption to modelize real noises. 
To answer this question, at least it would be necessary 
to calculate the first contributions of the finite value ex- 
ternal noise parameters. If this corrections are controlled 
then one can get confidence on the simplified assumption 
of white noise. 

Since the early work of Schldgl |l3j on the effects of 
fluctuations on a chemical interface, an intensive work 
has been devoted to describe the related problem of front 
propagation in the presence of an external noise source 



^UJfjjTjlPI . A complete study for the case of a white 
noise was presented in |^^] in which the front veloc- 
ity and its diffusive dispersion behavior was computed 
in terms of the effective white noise intensity. It was 
found that the velocity of the front increases with the 
noise intensity due to a systematic contribution to the 
kinetic terms. Actually the dispersion is subdiffusive for 
the so called pulled fronts, as recently has been shown 
[fijl. A variety of approaches has been adopted, from 
projection techniques P,fjJ to scaling arguments jl5), in- 
cluding a stochastic version of the multiple scale analysis 
as well as a non-systematic noise expansion 0] . The last 
one may be the most simplified one which grasps the sys- 
tematic contribution of the noise to the dynamics of the 
system. This contribution usually appears as a renormal- 
ization of the reaction parameters and it is the origin of 
the well known shift of the front velocity. 

The general aim of this paper is to find what are the 
most relevant effects of a real noise on two different mod- 
els which exhibit front propagation. We will see that for 
a fixed noise intensity, the noise correlation time is a rele- 
vant parameter which interpolates the results of the white 
noise limit with those of the deterministic case, but the 
role of the correlation length is different. 

Here, we will derive an analytical expression of the ef- 
fect of a spatio-temporal structured noise on an extended 
system governed by a Langevin reaction-diffusion equa- 
tion with multiplicative noise. We will closely follow the 
guidelines settled in for a one variable system, and 
also those in |l7j for multivariable system to deal with 
non white noises in extended systems. Our main differ- 
ence is that our analysis is done in the continuum space, 
and also that we present a more simplified way to get 
the first order contribution of the noise in the correlation 
time r which avoid the integration of a response function. 
Here we note that the continuum Langevin description 
may perfectly be adequate for describing reactive fronts 
even though chemical systems are discrete in nature, as 
was shown in Jl8[ . 

Our theoretical predictions have been applied to two 
systems: the Ginzburg-Landau and the Schlogl models. 
In the first case the noise induces the front by control- 
ling the stability of the new state versus the other unsta- 
ble steady state. In the second case, the noise does not 
change the steady states but controls its dynamics. 

The outline of the paper is as follows. Section II con- 
tains the main theoretical results and a discussion of some 
limiting cases. There we present the derivation of an ef- 
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fective dynamical equation which grasps the systematic 
contribution of the different noise parameters. In section 
III, we apply these results to obtain explicit predictions 
for the two models already mentioned and we discuss the 
numerical results obtained for these them and their com- 
parison with the analytical predictions. In section IV 
we summarize our conclusions. Several Appendices are 
devoted to technical aspects of our analytical methodol- 
ogy and the implementation of a particular algorithm to 
generate a patio-temporal structured noise. 



II. EFFECTIVE DYNAMICAL MODEL 

We consider the following stochastic partial differen- 
tial equation as a representative description of reaction- 
diffusion systems under multiplicative fluctuations: 

= c(^(x, t),d x ,a) + e l ' 2 gWMx, t). (1) 

where £ is a reaction-diffusion operator which explicitely 
reads, 



C(ip(x,t),d x ,a) = D 



d 2 ^{x,t) 
dx 2 



+ /GM)- (2) 



f(ip, a) and g(ip) are the reaction term and the coupling 
term with external fluctuations respectively, and 77(2, t) 
is a Gaussian spatio-temporal structured noise with the 
following statistical properties, 

(n(x,t))=o 

( V (x',t')r 1 (x ) t))=G(\x-x'\,\t-f\) = 

= C(\x-x'\) 7 (\t-t'\). (3) 

Also, for simplicity but not strictly necessary, we have 
made the assumption that this correlation function fac- 
torizes in a spatial and temporal part. To fix the nota- 
tion, and following the commonly accepted generic pre- 
scription Q , we define the three parameters of the noise, 
intensity, correlation time and correlation length, as fol- 
lows: 



ds / drG(r,s) 



■R 



r = — I ds dr G(r, s) s 
o- Jo Jn 



(4) 



Y = — J ds l dr G{r, s) r 2 



n 



We pursue here to find the systematic and most rele- 
vant effects of this type of noise. In general, the noise has 
two important effects, systematic and fluctuating ones, 
which cannot be exclusively associated with the deter- 
ministic and stochastic terms of Eq. (|l|) respectively. In 
fact noise acts in two different scales 0. Fast fluctu- 
ations in a short time scale modify the front shape and 



thus producing an effective front with different determin- 
istic properties. On the other hand, the slow fluctuations 
are responsible of the diffusive dispersion of the front po- 
sition . 

A naive way to get these systematic effects of the fluc- 
tuations is by analyzing the noise term in (Q) . Due to de 
multiplicative character of the noise, although 77 has zero 
mean, this is not the case for this stochastic term, 



(5) 



which, as a consequence, will give a net contribution to 
the dynamics. This can be explicitly shown in the fol- 
lowing way. By adding and substructing < I ) ({V'}) to our 
original dynamical equation (Q), we can write this equa- 
tion as 

= cMx, t), d x ,a)+ + e V2 W) X) t)t 

e 1 ' 2 K(il>, x, t) = eWgtyWx, t) - $({</>})• (6) 

This dynamics is statistically equivalent to the original 
one. Note that for the new noise term it is (H(ip, x, t)) = 
and it has a correlation which can be developed in 
powers of e 1 / 2 . We make now the Ansatz that if the noise 
allows for a well definite front structure, its systematic 
behavior will be described by the deterministic equation 



dip(x,t) 



= £W>0M),^,a) + $(M) 



(7) 



called the effective dynamics. In appendix A and B we 
present a detailed calculation of $ for small r, which is 
given by 



$(W) = eDC(0)r S ' 9 "( 
where 



dx 



h(^(x,t)) (8) 



h(iP(x,t)) = fWx,t),a) + 
e{C(0) + DC"(0)T}g(i>(x,t))g'(i>(x,t)) 
- e C(0)t 9 '(iP(x, t)) {gMx, t)), f(${x, t))} (9) 

is the new effective reaction term. The brackets are de- 
fined as 



{9 » / } = 9' f ~ 9 f ■ 



(10) 



and the primes on f(ip) and g(ip) indicate the derivative 
with respect to ip. 

Thus, we have ended up with one of the most impor- 
tant results of this paper, Eq. ([?]), which contains the sys- 
tematic contribution of the noise to our original dynamics 
(|l|) up to first order in r. In this paper will not study 
the effect of 7Z. This term is only relevant for those non 
systematic effects of the noise such as, for example, the 
dispersion of the front. The dependence on the parameter 
A is included in C(0) ~ <r 2 A -1 and C*"(0) - a 2 \~ 3 . As 
will be seen below, C(0) is the most relevant quantity. 
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Thus the main effect of the correlation length through 
C(0) is trivial. For this reason we will pay more atten- 
tion to the non-trivial influence of the term C"(0) fixing 
C(0) independent of A. 

A first check of the previous results will be provided 
by considering the better known case of temporal white 
noise in a lattice. Here we will first define the proper 
limit by which (^|) becomes a temporal white noise, and 
then see if ([?]) correctly reproduces the results found in 
Ref. §. 

The temporal white noise in time has a correlation, 



g(ip(x,t)) 



da 



(17) 



(r)(x', t')r)(x, t)) =2C(x- x') S(t - t') , 
where the spatial white noise limit is given by 
lim C(x-x') = S(x-x'). 

In a one dimensional lattice this takes the form of 



lim Ca = -P- 

A^o 1 Ax 



(11) 



(12) 



(13) 



In this limit one can see that all the integrals in (AS) 
vanish except the first one. The systematic dynamics of 
([?]) is then that of the effective reaction term given now 
by (a 2 = 1) 



MV-OM)) = f(ip(x,t),a) + e(0)g{i>(x,t))g'(i>(x,t)), 

(14) 

with e(0) = eC(0) = e/Ax. 

In this way we have recovered the results found in Ref. 
|^| . Note that one cannot consider right from the begin- 
ning a white noise in space because the ill defined 6(0). 

The value of C"(0) for a spatial white noise in the 
lattice is evaluated as, 



C"(0) 



C(l) -2C(0) + C(-1) 
(Ax) 2 



(Ax) 3 



(15) 



where C(±l) = has been used. For the case of a sp atial 
structured noise with A finite, all the integrals in ( A8) can 
be evaluated. 



III. APPLICATIONS AND NUMERICAL 
RESULTS 

We will now study the effects of a colored noise for 
two particular types of couplings g(<j>): a linear and a 
nonlinear one, which correspond to the Ginzburg-Landau 
and Schlogl models respectively. The noise will enter in 
the standard way Oil] as small fluctuations of the control 
parameter a, 

a -> a + e 1/2 i](x 1 t), (16) 
and thus the Langevin-type coupling function is given by, 



Numerical simulations of the Eq.(^) for the different 
models have been performed in a one-dimensional lattice 
of mesh-size Ax = .5. The length of the system is L — 
600. We have used a Heun algorithm with a time step 
At = .01. In all cases, D = 1 and a = —.1, except for 
the Schlogl model, where different values of a have been 
used. 

The noise is generated with a spatial and temporal 
structure as a Gaussian random numbers at each lattice 
point. The correlation function factorizes as in Eq. ([?]). 
The temporal part has a exponential decay (Ornstein- 
Uhlcnbeck process) with a correlation time r, while the 
spatial correlations have a triangular shape with a corre- 
lation length A. The numerical implementation of such a 
noise is described in appendix C. 

The initial condition for the Ginzburg-Landau model 
is a small pulse with a height of 0.01 and located at the 
middle of the spatial domain. In this way the initial per- 
turbation will spread off as two fronts propagating in op- 
posite directions. For the Schlogl model the initial field is 
a step-like function of value ^(0, 0) = 1 and ip(x,0) = 
in the rest of the spatial domain. The numerical calcu- 
lation of the mean front velocity and the steady state 
behind the front have been done as in Ref. EOjl. 



A. Linear coupling: The Ginzburg-Landau model 

This model have already been considered in the context 
of noise-induced fronts [^0| . We will now study how that 
picture is modified by a spatio-temporal structured noise. 
For this model the kinetic term is, 



/(</>, a) = -V(o + V 2 ) 



(18) 



and as a consequence, the noise coupling term is linear, 
9{<t>) = 0- (19) 

We will have then that the effective dynamics given by 
Eq. © is 

^ = D d -^ + Hi>(x,t)) (20) 

with a new kinetic term ([)]) given now as, 

h(ijj(x,t)) = -i/>(a' + &> 2 ), (21) 

where the new kinetic parameters are, 

a' = a - e{C(0) + DC"{0)r} 

b' = 1 + 2eC(0)r. (22) 

Following the linear marginal stability criteria J21[ , the 
velocity of a pulled front is controlled by the linear term 

as, 
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= 2jD(-a + £ (0) (1 + ^Spr)). (23) 

Note that a' < in order to have a front. This result, 
however, has been deduced for a small enough r. Never- 
theless, we can conjecture a generalization of Eq.(|23|) for 
any value of r considering that the values of the velocity 
for r — (temporal white noise limit) and for t = oo 
(deterministic case) are known. In this way, the simplest 
regularization of Eq. (^3|) , which is a monotonous function 
on r, is 



1 



D(-a 



e(0) 



DC"(0) 
C(0) 



(24) 



Moreover, Eq.(J24|) generalizes the temporal white noise 
result in terms on a renormalized noise intensity defined 
as, 



6(0) 



DC"(0) 



(25) 



which does not present any singularity because always 
C"(0) < 0. 

Taking e(0) = eC(0) as a constant we have the follow- 
ing behavior. For increasing r (A fixed), e-n decreases, 
and we arrive up to the deterministic value of v. Never- 
theless for fixed r and increasing A, C"(0)r/C(0) ~ t/A 2 , 
and then c-r, increases, and as a consequence, v increases 
up to the temporal white noise limit. This is a non trivial 
effect of A which needs a finite value of r to appear. 

As already discussed in [po| , another important quan- 
tity in this model is the field behind the front which is 
induced by the noise, and thus it is highly fluctuating. 
From Eq.(pl|), the homogeneous deterministic stationary 
value behind the front can be calculated as, 



(26) 



As we don't know how the higher order corrections on 
b' are, we expect a poorer agreement for ip st than for the 
velocity. However, we can get an idea of the relevance of 
V by numerically inspecting the quotient of v* over -0 s t 
which depends on b'. Indeed, from Eqs. ( |23| ) and (26), 
this dependence is 



2 V, 



(27) 



The analytical predictions (|24|), ( |26|) and (27) are im- 
portant results of this paper that will be checked numer- 
ically. Due to the different role of the noise parameters 
t and A we will discuss two cases separately. 



1. Spatial white noise in the lattice 

For this case all the simulations perfectly agree with 
the theoretical results of Eq. (p4[). In fact, fixing the 
noise intensity e(0) and increasing r the mean velocity of 
the front drops monotonously to the deterministic value 
(see Fig. @. All fi gures are in dimcnsionlcss units). 

Our analytical calculation (dashed lines) only describe 
the corrections to the white noise case at order 0(t). 
However, this can be considered quite relevant as the de- 
pendence of v and ip s t on r drops down very rapidly near 
t = 0, and our first order approximation succeeds to 
grasps this pronounced slope (see inset of fig.0). More- 
over, the extended analytical prediction (|24|) shows a very 
good agreement with numerical data for all values of r. 

With respect to the average mean field, the agreement 
is more qualitative (see Fig.0). The numerical results for 
the effective parameter b', evaluated from Eq. (|27]), are 
represented in Fig.^. They support the initial growth of 
b' predicted by the theory For greater values of the cor- 
relation time b' stays bounded by its deterministic value 
b' d = 1. This fact may explain why numerical values of 
ipst seem to depend only on the linear coefficient a'. 

From previous work [|20f| , we already know that ijj s t has 
a systematic error that slightly increases with noise in- 
tensity. Hence, the deviation found here (Fig.||) is not 
due to the presence of temporal correlations, but a prob- 
lem related to the fact that we are measuring a highly 
fluctuating quantity as it is ip s t- In any case, our theoret- 
ical prediction is consistent qualitatively with numerical 
simulation results. 



2. Spatio-temporal structured noise 

To study the non trivial effects of a finite correla- 
tion length on the dynamics of the front, we have to 
pay attention to the effects coming from the quantity 
C"(0) ~ a 2 /\- z . 

In Fig.(|4|) we can see the front velocity versus r for 
different values of the correlation length A of the noise. 
Continuous lines correspond to the analytical prediction 
(pi|). As can be seen, for a finite correlation length the 
agreement is only qualitative, and improves for noises 
not too much away from a spatial white noise (triangles) . 
This behavior can better be appreciated in Fig. (||) where 
we have plotted the velocity versus A for different values 
of the correlation time. 

According to our definition of A (CS), we have that 



A = for the spatial white noise in the lattice. Here it 
can be seen that our analytical scheme may qualitatively 
describe the effects of a finite A only for small values of 
r and A, which is not the case for A = where the accor- 
dance is very good. For a finite A, we have also observed 
a clear departure from the analytical results in the case of 
the dependence on r of the mean stationary value of the 
field behind the front, as well as for the cubic coefficient. 
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Also, the numerical results show a systematic decrease 
of the velocity when increasing the correlation length for 
forr = 0, which is not predicted by our analytical results 

On the other hand, for t > 0, the velocity tends to 
grow with A for small correlation lengths as the theory 
predicts. Indeed, this agrees with what we have observed 
in preliminary numerical simulations for a quenched 
white noise. Hence, for increasing A (r fixed) the nu- 
merical results suggests a non monotonous behavior of 
the velocity, which may increase at small A, but would 
always decrease at long correlation lengths. 

We have not found yet an explanation for this effect. 
We believe that there is a interplay between the corre- 
lation length of the noise and the typical length of the 
front, which is given by its width. Indeed, this effect 
could be related to the observed distorsion of the leading 
edge of the front and the possible formation of a prefront 
in the presence of a large spatial correlation length of 
the noise. Then our initial assumption of a well defined 
mean front profile is not fulfilled and, as a consequence, 
the theoretical scheme cannot be applied. 

For the sake of completness, in this figure we also show 
what are the trivial effects of of a finite correlation length. 
In this case, the noise has been generated such that the 
noise intensity a 2 , defined in (^), remains constant. This 
can be acomplished by changing the previous whcighting 
function g$ by a factor (V2m + (See appendix C). 
Squares correspond to the numerical results (L — 2400). 
The dashed line is the theoretical prediction ( pi| ) for 
this case. As can be seen, there is a monotonous de- 
cay of the velocity because now the most dominant term 
e(0) ~ <7 2 /A also decays with A. Thus, the qualitative 
behavior is completely different. This confirms that the 
previous studied dependence on A correspond indeed to a 
non trivial effect of the correlation length. Note that our 
theoretical scheme succeeds better in describing quanti- 
tatively the trivial effects of a finite A. 

B. Nonlinear coupling: The Schlogl model 

The general model was introduced by Schlogl in C3] 
when studying the fluctuations of an interface. Here we 
will consider a particular version of it that was studied 
in the presence of an external white noise in Ref. It 
corresponds to the reaction term, 

/(^,a) = -rli(j> + a)(^ - 1), (28) 

which implies a nonlinear coupling with the noise, 

g{il>{x,t)) = VW> - !)• (29) 

Taking into account these definitions, the effective de- 
terministic part of Eq.(^|) becomes, 

dip(x,t) d 2 ip(x,t) ,,,,, 



+ Z )e( 0)r2 ( 2^-1)(^) 2 (30) 
where the reaction term is, 

h(ip(x,t)) = a'ip + b'ip 2 + c'ip 3 + 5d> 4 - 2d> 5 , 

(31) 

with the effective kinetic parameters, 

a' = a + e-jz 

b' = 1 - a + d! - 3 en 

c' = -1 - 4d' + 2e n 

d! = e(0) r, (32) 

where e-jz was defined in Eq.([25|). 

Our point of interest in this model is the mean front 
velocity because the steady states for the front ip = 0, 1 
are not modified by the noise. Due to the prefactor of 
the KPZ-like term, for tp > 1/2, any deviation from the 
homogeneous state ip = 1 tends to grow, while this is op- 
posite at points where ip < 1/2 for any deviation from the 
state ■0 = 0. Thus, the effect of this term is to shorten the 
width of the front, i.e., to select a greater decay mode of 
the front, thus slowing down its propagation pljj . Hence, 
the expected slowing down of the front due to temporal 
correlations of the noise arises here in two ways: The 
first one is by means of the usual renormalization of the 
coefficients of the reaction term f(ip). The second type 
of corrections come from the new KPZ-like term. 

We expect that this front will exhibit the general 
regimes of front propagation |2l|| , i.e., a linear (pulled 
front), a non- linear and a metastablc (both, pushed front) 
regimes. In the linear regime, the velocity depends only 
on the linear coefficient of the reaction term a', and it is 
then given by Eq. (pi|). 

While a crossover from a metastable to a nonlinear 
regime is trivial to determine, being nothing more than 
condition a' = 0, the transition between the linear and 
the nonlinear regimes is far more complicated to locate. 
This calculation of this point requires an analysis of (all) 
higher power terms of the reaction, determining a com- 
plete solution in the comoving system and requiring then 
that the asymptotic behavior is such that the coefficient 
of the slowest decay mode vanishes. 

For the case of a temporal white noise, as was found 
in HU, there is only a renormalization of the parame- 
ters of the kinetic terms, in such a way that the effective 
dynamics is equivalent to the deterministic one up to a 
rescaling of the coefficients. Hence, the location of the 
different regimes can be directly determined from those of 
the deterministic case. Unfortunately, to our knowledge, 
the first procedure is hopeless for Eq. (|30|) . Neither can 
this dynamic directly be compared with the determinis- 
tic case. However, the transition between linear and non 
linear regime is always continuous and, as we can cor- 
rectly describe the linear regime, this fact will help us to 
numerically locate the transition for this model. 
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Nevertheless, there is still some hope for an analytical 
prediction. The type of dynamics given by Eqs. (|30|)- 
( |32| ) usually are relevant near the transition point a' = 0, 
where the dynamics given by Eq. (0) can be simplified by 
means of an amplitude expansion. In this case, and as 
long as the noise intensity is low enough, our effective 
equation would also lie near threshold (a' ~ 0). Assum- 
ing this situation, the spatial variations of the field take 
place on a typical length scale of order q^ 1 = ■sjDja. 
A crossover between non linear and linear regime means 
that the nonlinear terms start to dominate the growth 
rate off the initial steady state. Thus, this transition 
takes place when both, a! ip and b' ip 2 are of the same 
order of magnitude. This will be the case for 

ip ~ b' ; b' 2 ~ a'. (33) 

Then the KPZ and the ip 4 terms, both will be of order 
ffp, while the term ip ( ^^'^ J will be of order . Hence, 
near threshold, only the first three terms of the effective 
reaction ( [31] ) will be relevant. But this equation is just 
the standard Schlogl model that is exactly solvable. This 
will have a sense only if it is d < .which we will assume 
to be so due to its expression in (p2|). As stated above, 
we will also require b' > 0. The new stationary states of 
this approximation are given by 

For a' d < and b' > 0, it is ip+ > ip-. If we write then 
the reaction term as 

h cr W(x, t)) =diP{il>- 1>+) & - (35) 

the velocity of a front connecting ip = ip + and ip = is 
given by Q 

v nl = ^2dD{t± - ^_}. (36) 

Given the results found for the previous model, we 
expect that our system (^Tj) will present fronts in the 
linear regime for a high enough noise intensity e(0) and 
small enough r. 

Although we do not know the corrections of d! beyond 
O(r), we expect the velocity to start diverging from 
for some finite value of r. But, by adequately choosing 
the parameters, the transition between the linear and 
nonlinear regimes can be obtained at small enough val- 
ues of t for our approximations to be applicable. Hence, 
for values of the a 1 near threshold, and within a neigh- 
borhood of the crossing point, the nonlinear velocity v n i 
will approximately be given by Eq.([36"|). 

We will show that numerical simulations support this 
analysis. In Fig.(^) we have plotted the numerical results 
of the front mean velocity versus r. For all three plots, 
the value of the deterministic linear coefficient a is such 



that the deterministic front lies right within the nonlin- 
ear regime. While for the white noise case all fronts move 
within the linear regime, the one for e(0) = .1 is only 
marginally inside. Increasing then r, the effective linear 
coefficient a' decreases. For small values of the correla- 
tion time, the fronts will still lie within the linear regime, 
except for the e(0) = .1 case, for which the front enters 
the nonlinear regime immediately for any finite value of 

T. 

Thus for small values of r the front moves with the 
linear velocity. We can see that our first order approx- 
imation (dashed lines) also reproduces for this system 
the initial steep fall of v. Our analytical continuation v* 
shows up a perfect agreement with the numerical results. 
By further increasing the correlation time we can shift 
the front into the non-linear regime. 

For this one we only have a rough approximation for 
the front mean velocity given by d3q) (dot-dashed lines). 
This approximation is valid only near the critical point 
a' = and gets worse as we move away from it. Although 
our analysis gives us only the corrections due to a finite 
r up to O(t), with the renormalized noise intensity (|25|), 
and staying close to the critical point, we obtain a sur- 
prisingly good analytical estimation for the non linear 
velocity up to values of r of order 0(1). One can no- 
tice that our predictions are just O(r) seeing that the 
velocity diverges from the deterministic value at r high 
enough. This is to be expected as in obtaining (^) we 
have used only an O(r) approximation for the squared 
b'tp 2 and the stabilizing term di/j 3 , with both, b' and d 
linearly increasing with r for high values of r. But this 
dependence is obviously incorrect as a infinite energy dif- 
ference between two metastable states would give rise to 
an unbounded propagation velocity, which is not the case. 

The estimation of the crossing point between the linear 
and the non linear regimes, given by the crossing point 
between the two curves v* nl and v* , is consistent with the 
numerical results. At that point, these show up a pro- 
nounced and increasing departure off the linear velocity 
v*. For the e(0) = .1 case, the departure from the the- 
oretical result v* starts already for r > 0. This can be 
best seen in the inset of figure ^]. 

IV. CONCLUSIONS 

We have studied a general reaction-diffusion system 
which exhibits fronts in the presence of spatio-temporal 
structured external noise. We have derived an effective 
deterministic dynamical equation for the front which con- 
tains the main effects of the noise. These show up as 
renormalization of the original reaction terms of the de- 
terministic system and also new terms. The relevance 
of those effects are discussed for two prototype models: 
The Ginzburg-Landau model (noise linear coupling) and 
the Schloegl model (noise nonlinear coupling) . Although, 
our analysis is valid only for small correlation times, we 



G 



are able to grasp at least the main features of the effects 
of a time colored noise on an extended system, namely 
that there is a pronounced slow down of the front velocity 
already for small but finite correlation times r. 

We have obtained an analytical expression for the front 
velocity in the linear regime that we have extended to be 
valid for any value of the correlation time r of the noise. 
The numerical simulations present an excellent agree- 
ment with the analytical results in this regime. For the 
non linear regime, we can only give an approximate ex- 
pression for the velocity, when the system is near thresh- 
old. Although it still is valid only for r not too large, it 
goes beyond the first order approximation considered in 
our analysis when the system is brought close enough to 
threshold. 

Finally, we have obtained the non trivial influence of 
the noise correlation length. Our results here are of a 
more limited validity, and this case needs a further study. 
Nevertheless, our numerical results suggest a non triv- 
ial behavior of the velocity when varying the correlation 
length of the noise, with an increase of the velocity with 
A for finite r and small A. This would be a novel result, 
because previous reported studies evaluate only the triv- 
ial dependence on A of the front velocity, slowing down 
with increasing correlation length |2^] . 

Hence, we have shown that our procedure of separating 
the systematic contribution of the noise from the origi- 
nal dynamics gives reliable information for front dynam- 
ics in the presence of spatio-temporal structured noises. 
The systematization of the present procedure, if possible 
at all, and its application to other situations would be 
extremely interesting. 
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APPENDIX A: ANALYTICAL DERIVATION OF 
THE SYSTEMATIC EFFECTS OF A 
STRUCTURED NOISE 

The systematic contribution (<f>(x,t)) of the noise is 
given by 



(Al) 



This average con be calculated by using Novikov's theo- 
rem in the following form, 



=e 1/2 / dx' dt'G(\x-x'\,\t' -t\) 
Ju Jo 

Sip(x, t) 



g'(il>(x,t)) 



5r)(x',t') 



(A2) 



Hence, the determination of ($>(x, t)) reduces to that of 
the response function 



Q(x,x';t,t') 



Sip(x, t) 
6r)(x', V) ' 



(A3) 



Following |l6|] and jl7J , we will consider the contribution 
of the noise at first order in r in the approximation of 
small r. This means that temporal correlation decays 
very strongly for t' ^ t leaving relevant in QB3| ) only the 
values of the integrant for t' close to t. Thus we may 
expand Q(x,x';t,t') in powers of (t' — t) around t = t 
and take all up to the first order 



Q(x,x';t,t') = Q(x,x';t, t) + 



dQ 
dt 1 



(f - t) + 



(A4) 



In appendix B we present a detailed derivation of this 
second term. 



Now we can rearrange (A2) in two terms, the first one 
being the zero order or white noise contribution, i.e, the 
one we get in the limit r — > for fixed e and A, while 
the second one represents the contribution of the colored 
noise at first order in r, 



(*(x,t)) = (*oOM)> + ($iOM)>- 



(A5) 



Recollecting relations (A2),(AJ) and (B4),(B7), and after- 
calculating the spatial integral, we obtain 

($„(*, t)) = e (gty(x, *)) g'Wx, t))\ / dt' G(0, \t-t'\) 



($ 1 (x,t))=e(g'(ilj(x > t)) 

fMx,t))( d ^ x >^ 2 



-Dg" 



dx 



{g{i>{x,t))J{i>{x,t))} (A6) 



dt'G(0, \t-t'\) (f-t) 



Dg(ij(x,t)) / dt'G"(0,\t-t'\)(t'-t) 



where the primes on G(0, (t — t')) indicate derivatives 
of G((x — x'), (t — t')) with respect to x' , evaluated at 
x' = x, and 

{g(iP(x,t))J(iP(x,t))} = g'(ip(x,t))f(ip(x,t),a) - 

- g(il>(x,t))f'Mx,t),a). (A7) 

As we are interested in the approximation of small r, 
which amounts to consider observation times much more 
greater than the characteristic correlation time of the 
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noise, we can extend then the limits of these integral 
up to oo. 

At this point, further assumptions on the correlation 
function must be done in order to obtain any analytical 
prediction. Assuming that G(cc, s) factorizes as in (||), 
the above integrals can be written as 



a formal solution as it was in M] and |T^] for non spa- 
tially dependent (O-dimensional) systems. The term (B4) 
gives the contribution JAS| ) to the systematic effect of the 
noise. Expanding Q(x, x'; t, t') in powers of (t'—t) around 
t' = t and taking all up to the first order 



dt'G(0,\t-t'\) = (7(0) ~ -- 

A 

/ dt'G{Q,\t-t'\){t' -t) = -tC(0) ~ 4^ 
Jo ^ 

[ dt'G"(0,\t-t'\)(t' -t) = -tC"(0)~^. (A8) 

Jo * 

where the temporal part j(s) is considered to be normal- 
ized to 1. 

Finally, we can write 



Q(x,x';t,t') = Q(x,x';t,t) + 



(f-t) + 



(B5) 



The second term can be obtained by directly deriving 
(B3) with respect to t' . This gives 



dQ(x,x';t, t') 
dt> 



dsl(x,s)-{^^-,Q(x,x';t,t')} 

(B6) 



($>o(x,t)) = eC(0) (g(i>(x,t)) g'(iP(x,t))j (A9) 
(§i(x, t))=-e C(0) r ( g'(i>(x, t)) [{<?(</>(*, t)), f(^(x, t))} 



As long as we are interested in the limit t' — > t, the 
details of T(x, s) are not important for the last term in 
( p36|) vanishes in that limit since it is a regular function 
in s. Fort t' = t, and substituting equation ([!]), 



dx 

+ eDC"(0)r(g'(ilj(x,t))g(ij(x,t)) 



(A10) 



dQ(x,x';t, t') 



Of 



dQ(x,x';t',t') 



APPENDIX B: RESPONSE FUNCTION 

The determination of (&(x, t)) reduces to that of the 
response function 



dip(x,t') 

C^(x,t),d x ,a) + e^^M))^,*)^ - 
fd£, 



{x, t) + eWg'Mx, t)) r](x, t) Q(x, x';t, t) . (B7) 



Q(x,x';t,t')- 



(Bl) 



6r]{x',t') 
The meaning of (|l]) is given by 

/t+At 
ds C(tp(x, s), d x , a) 

rt+At 

e 1 / 2 / dsg(ip(x,s))r,(x,s), (B2) 



and as long as r > 0, the last integral is well defined as 
a Ri emann integral. By taking the functional derivate of 
(B2) with respect to the noise rj(x',t') we get 



Considering the initial condition (B4) and substituting 
the expression of the nonlinear (differential) operator C 
given by (El), the last relation reduces to 



dQ(x,x';t,t') 



Of 



g'(il>(x,t))f(ilj(x,t),a) 



g(i>(x,t)) f'(ij;(x,t),a)\6(x' - x) + 

De^{g'Wx,t))(^^)S(x'-x) 
d 2 



dx 2 



g(iP(x,t))5(x' -x) }. (B8) 



Q(x,x';t,t') = Q(x,x';t',t') + 
J ds{ — {x,s) + e 1/2 g'(iP(x,s))r ] (x 1 s)}Q(x 1 x , ;s,t') 
; t > t' 



The second order derivative of the last term gives 

d 2 i>(x,t)^ 



g"Wx,t))(^^l) 2 + g'(^x,t))C- 



dx 2 



with 



Q{x, x'; t, t) = e 1/2 g(ip(x, t)) S(x' - x). 



(B3) 



(B4) 



S(x -x) + 2g {tp(x,t)) ( dx ) ^ 

d 2 8(x' - x) 



-g(tp{x,t)) 



dx 2 



(B9) 



Equation ( JF33| ) is an integro-differential equation for the 
response function for which it has not yet been found 



With this result in mind, the Laplacian terms in (BS) 
mutually cancel, while the terms proportional to the S 
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give rise to the first integral in ( A10). Taken into account 
the following relation of the derivative of a S 



_0_ 

Oa- 



Six' - x) 



d_ 

dx 



7 8(x' - x) 



(BIO) 



the contribution to (|A2|) of the first and second order 
derivative of the 5 in (B9) will be 



dt' 



2G'(0,(t-t'))g'^(x,t)) 



dip(x, t) 
dx 



G"{Q,{t-t'))g(i>{x,t)) (t'-t). (Bll) 



This result give rise to the last terms in (A10), where 
the contribution proportional to G'(0, it — t')) has been 
discarded because of the spatial isotropy of the noise. 



APPENDIX C: NUMERICAL ALGORITHM FOR 
GENERATING A SPATIO-TEMPORAL 
COLORED NOISE 

Here, we will define a spatio-temporal structured noise 
that is very simple to implement numerically, and which 
is the one we have used in this work. This type of noise is 
obtained by rewriting the spectral method M as a linear 
transformation of a more simple noise field in real space. 

We define our spatio-temporal colored noise in each 
lattice cell i and at time t, as 



7?i(i) = Ax Vj(t)9 



(CI) 



Being G already factorized, the linear transformation 
( |Cl|) assures that the correlation function of %(t) will 
be of the desired form (||). As g is arbitrary, we are free 
to impose the condition that the value of Go(0) equals 
that of Go(0), i.e, 



G o (0) = G o (0) 



(C6) 



Here we are interested in gi having a finite range, for 
simplicity, we assume that <?.; is a constant g inside the 
interval — m < i < to, but zero otherwise. Then the 
condition (C6) implies gi — ( Ax y^2m + l) )^ 1 .Now 
it is a simple calculation to show that rji{t) is a spatio- 
temporal structured noise with a correlation, 



Gi = C, 7(|*|) 



G 



(2 m + 1) 



[2 m 



1 



0(2m-|Z|) 7 (|*|). (C7) 



At equal lattice points this function decays exponen- 
tially in time, and at equal times it has a triangular decay 
as a function of lattice point difference. From this ana- 
lytical expression it is straightforward to obtain the noise 
intensity and correlation length, 



a 2 = (2m + l)a 2 v 



A = \j — m (m + 1) Ax, 



(C8) 



being the correlation time r. One can check now that in 
the lattice white noise limit to = and a 2 = er^, and 
then A = 0. For to — 1, we get that A = 1.15... Ax. 



where the index j labels a domain of cells around the cell 
i, and gi is a weighting distribution with the isotropic 
property 

.9-* = 9i- (C2) 

fjiit) is an Ornstein-Uhlenbeck process in the lattice 
cell i, statistically independent of the other lattice points 
(white noise in space) . Its value is generated through the 
linear Langevin equation, 



at t t 



(C3) 



in terms of a Gaussian white noise with a correlation, 



Ut%it)) = 2*l-£-5it-t'). 



(C4) 



In this case the correlation of the noise to(£), is given by 
(fj l+j is)fj j iO))=G l is) = C ll i\s\). 



7(N) = -f e" 
= Ax' 



(C5) 
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FIG. 2. Mean stationary value of the field versus t. The 
inset is a amplification of the small r domain. See previous 
figure for the symbol notation. 



1.2 



1.0 " 



1 L-_ 



1 2 3 4 5 

FIG. 3. Cubic term coefficient of the Ginzburg- Landau 

model as the quotient of v to <f> versus t. Notation is the 
same as in previous figures. 



FIG. 1. Front mean velocity versus noise corelation time 
t for the Ginzburg- Landau model in the presence of a Orn- 
stein-Uhlenbeck noise in time, white in space. Values of 
the parameters: Circles correspond to the deterministic case, 
v* im = .62; Triangles, to e(0) = .2, and squares, to e(0) = .05: 
dashed lines are the analytical prediction up to 0(t) (p^|), 
whereas continuous lines represent the corrected prediction 
See text for the values of the other parameters. 
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2 3 
x 

FIG. 4. Front mean velocity versus r for different correla- 
tion lengths. Here it is e(0) = .2. 



v(x,X)o.9 - 




v(x) 1.4 




FIG. 6. Front mean velocity versus noise corelation time r 
for the Schlogl model. Circles correspond to the deterministic 
case, sim = 1.131; dashed lines are the analytical results 
up to 0(t), whereas continuous lines represent the analytical 
corrected results v* Eq. (jj^). Dot-dashed lines correspond to 
the prediction (p^j. The inset amplifies the domain of small 
t . By tuning t, the front shifts from a pulled regime to a 
pushed one. This crossover corresponds to the points first 
leaving the theoretical curve v* and are approximately given 
by the arrows. 



FIG. 5. Front mean velocity vs. A for different correla- 
tion times. e(0) = .2 for hollow triangles. Filled up-triangles 
correspond to the white noise case, while circles to the deter- 
ministic case. Squares correspond to the trivial noise effects 
with fixed noise intensity for r = .05 and the long-dashed line 
is its theoretical prediction (see text). Each set of verti- 
cal points along the horizontal axes correspond, in increasing 
order of X, to m — 0,1,2,3,4, respectively. See appendix C. 
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